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1 . INTRODUCTION 

To be successful in the development of 
satellite retrieval methodologies for the 
determination of cirrus cloud properties, we must 
have fundamental scattering and absorption data on 
nonspherical ice crystals that are found in cirrus 
clouds. Recent aircraft observations (Platt et al . 
1989) reveal that there is a large amount of small 
ice particles, on the order of 10 pm, in cirrus 
clouds. Thus it is important to explore the 
potential differences in the scattering and 
absorption properties of ice crystals with respect 
to their sizes and shapes. In this study the 
effects of nonspherical small ice crystals on the 
infrared radiative properties of cirrus clouds are 
investigated using light scattering properties of 
spheroidal particles. 

In Section 2, using the anomalous 
diffraction theory for spheres and results from the 
exact spheroid scattering program, efficient 
parame terizatat ion equations are developed for 
calculations of the scattering and absorption 
properties for small ice crystals. 
Parameterization formulas are also developed for 
large ice crystals using results computed from the 
geometric ray- tracing technique and the Fraunhofer 
diffraction theory for spheroids and hexagonal 
crystals. This is presented in Section 3. 
Finally, applications to the satellite remote 
sensing are described in Section 4. 


2. PARAMETERIZATION FOR SMALL ICE CRYSTALS 

Because the effects of nonsphericity would 
be small in the infrared region, small ice crystals 
may be approximated as spheroidal particles. For 
small size parameters of spheroidal particles, 


scattering parameters computed from the Mie theory 
are adjusted to match with those computed exactly 
by the method of Asano and Sato (1980) based on the 
anomalous diffraction theory (Latimer 1980) . 

For a spheroid of semi -axes a t and va fl the 
equivalent spherical radius is given by 

4g * a e [ v/(cos 2 ( ♦ v 2 sin 2 C) 1/2 ) kl . (1) 

where { is the angle between the incident direction 
and the rotation axis of the spheroid, a e the 
radius of the revolved circle, and v the axial 
ratio. The extinction or scattering cross section 
of the spheroid is given by 

C<{) - <?.ph «**[(cos 2 ( ♦ v 2 sin 2 C ) 1/2 ] k2 , ( 2 ) 

where Q ftph is the optical efficiency factor from 
the Mie theory for a sphere of radius, a g , and 
refractive index, m (- m r - 1 m^). If k^ is 1 in 
Eq . (2), the right-hand side becomes Q gph times the 
geometric shadow area of the spheroid. The cross 
section for randomly oriented spheroids is expres- 
sed by 

c..* - ll' Z c *.* (Osin c dc ■ (3) 

The unknown constants, k^ and k 2 , in Eqs . (1) and 

(2) are determined by matching C e /A, where A is the 
average geometric cross section of the spheroid, 
with that computed exactly by Asano and Sato (1980) 
as follows: 

*< - 1.0, Jc 2 « 0.98, for vzl 

l (4) 

■ 0.96, *2 * 1.08, forv<l.\ 


PARAMETERIZATION FOR LARGE ICE CRYSTALS 


Further, by comparing with values computed from the 
method of Asano and Sato, the extinction cross 
section for randomly oriented spheroids may be 
expressed by 


C % ♦ for <0.07 t 

C t , for JDi *0.07 . 


(5) 


where C # yeg is the absorption cross section for the 
volume equivalent sphere. The single scattering 
albedo is given by 


<«> 




C % * C Mf 


( 6 ) 


When the spheroid shape becomes flat or elongated, 
scattered light rays are inclined toward forward 
directions. This effect on the asymmetry factor is 
mimicked by adjusting the refractive index and 
size . 


With this procedure, the single scattering 
properties of randomly oriented spheroids can be 
efficiently calculated. Figures 1 (a) and 1 (b) 
show the single scattering properties of randomly 
oriented oblate and prolate spheroids with the 
axial ratio of 2 as a function of x G - 2*r G /X for 
11 and 12 pm wavelengths, where r G is the radius of 
area-equivalent spheres. The solid and dash-dotted 
lines are calculated by means of the preceding 
parameterization. Circles and crosses are values 
computed exactly from the method of Asano and Sato. 
The accuracy of the present parameterization is 
within about 3%. 


3. 


For large size parameters we use the 
geometric optics approximation to develop para- 
meterization equations. The single scattering 
properties of spheroids are computed from the ray- 
tracing method and Fraunhofer diffraction theory. 
Using these computed values , parameterized formulas 
for the single scattering properties of large ice 
crystals are derived. 

Figure 2 shows the scattering phase 
functions for randomly oriented oblate spheroids 
with the axial ratio, 1/v - 2, and m - (1.33, 0) 
computed from the ray-optics approximation and the 
exact method. The size parameter used is 30, which 
is twice the value used in Fig. 9 of Asano and Sato 
(1980). The peaks at the 0* scattering angle that 
are produced by the Fraunhofer diffracted light 
have about the same values for the two methods. 
The peak appeared at 0 - 76*, computed from the 
ray-optics approximation, corresponds to the 
maximum created by one internal reflection (primary 
rainbow feature). Since the rainbow angle is 
dependent on the incident direction on the 
spheroid, the rainbow peak Is smaller than that of 
spheres. Based on the exact method the rainbow 
angle is located at about 6 - 94*. The peak is 
smaller and broader and shifts toward backward 
direction. This behavior is the same as that in 
the case of spheres (Liou and Hansen 1971) and 
circular cylinders (Takano and Tanaka 1980) . 
Glory-like peaks appear at the backward direction 
(0 - 180*) on both methods. If the size parameter 
is larger, these two methods give closer values. 
This comparison verifies the ray- tracing program 
developed in the present research and the 
Fraunhofer diffraction approach developed by Takano 
and Asano (1983). Hence, it is possible to para- 
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Fig. 1: Comparison of single scattering properties for randomly oriented (a) oblate and (b) prolate 
spheroids between the parameterization (solid and dash-dotted lines) and the exact computation 
(circles and crosses ). The complex refractive indexes of ice used are m - 1 .0925 - iO .2480 at X - 
11 pm and m - 1 . 280 - 10 .4133 at X - 12 pm. 






diffraction 

formula for spheres , 

the 

single 

scattering 

albedo for intensively 

absorbing 

spheroids may be written in the form 




« - 1 - (1 - oj 2 Jt(cz) . 


(10) 

where 

*<*> - 4 * • 

w c 


(ID 

As a result, 

w - 1 when z - 0, and u 

- Q 

when z 

— *■ *>. To be 

consistent with Eq. (9), 

z is 

chosen 


as the independent variable. The constant, c, in 
Eq . (10) should be determined by fitting this 

formula with ray-optics computations. u in Eq . 
(10) is given by Q SCfl */2. Q sca «, is 1 plus the 
external reflected energy. Figure 3 (b) shows the 
single scattering co-albedo as a function o f z (> 
0.4). The constant in Eq . (10) depends on the 

axial ratio of spheroid. 


The extinction efficiency is given by 
<?.(*) - 2 * ^ (Q e (x c ) - 2) . 


( 12 ) 


Fig. 2: Comparison of the phase function for 
randomly oriented oblate spheroids using the ray- 
optics approximation (solid line) and the exact 
theory (dash-dotted line). The axial ratio, the 
complex refractive index, and the size parameter 
used are 2 f (1.33, 0) and 30, respectively. 


meterize the single scattering albedo and asymmetry 
factor of large spheroids using these programs. 

When absorption in a spheroidal particle is 
weak, the absorption cross section can be expressed 


by 


4TUDj 


(7) 


where m. is the imaginary index of refraction and 
r y the radius of the volume equivalent sphere. 
Thus, the absorption efficiency is given by 


«.* A 

* r B 


Aiun^ r y 


( 8 ) 


The single scattering co-albedo may then be expres- 
sed by 


1 - 0 ) - 


Q* 

T 


n*l 


* 


(9) 


Using values of o computed by the ray optics ap- 
proximation for a number of z, the coefficients, 
a n , are determined. Figure 3 (a) shows the single 
scattering co-albedo for randomly oriented oblate 
spheroids as a function of z. The solid line is 
calculated from Eq . (9), whereas symbols (o, + , and 
x) are computed from the ray-optics approximation. 
The fitting is sufficiently accurate in the range 
of z - 0.4. On the other hand, when absorption 
dominates, the single scattering co-albedo should 
approach to a constant value, which cannot be 
expressed by Eq . (9). Using the anomalous 


With this functional form, Q e - 2 when x « and 
Q e - Q e (x c ) at the crossover point, x £ . The cross- 
over point is the size parameter which separates 
small and large particle regions. The crossover 
point is determined in the following manner. 




Fig. 3: Single scattering co-albedo for randomly 
oriented oblate spheroids as a function of z for 
(a) weak and (b) strong absorption. The assumed 
axial ratios are 2,3, and 5, which are represented 
by o, and x, respectively . 






Figures 4 (a) and 4 (b) show the composite single 
scattering albedo for randomly oriented oblate and 
prolate spheroids with the axial ratio of 2 as a 
function of the size parameter, x , at a wavelength 
of 11 pm. The solid lines are the large particle 
parameterization using Eqs . (9) and (10), whereas 
the dot-dashed lines are the small particle 
parameterization using Eq. ( 6 ). For comparison 
purposes, the values computed by the exact method 
are represented by circles. From these figures the 
crossover point, x c , is about 22 in units of the 
area equivalent size parameter. A similar 
procedure can be applied to the asymmetry factor. 

For size parameters larger than about 30, 
parameterization equations have been derived using 
the results computed from the geometric ray- tracing 
technique for hexagonal crystals (Takano and Liou 
1989; Liou et al . 1990). The extinction cross 
section for randomly oriented hexagonal crystals of 
the same size can be expressed by 


c. 


=<")’ ijr. »> 


(13) 


where w is the width of the basal plane of a hexag- 
onal crystal and D the crystal maximum length. 
Based on numerical calculations involving a number 
of wavelengths, particle sizes, D, and aspect ra- 
tio, D/W, we find that the single scattering albedo 
can be parameterized in the form: 


1 - E . 

17*1 


for z < 0.4 , 
(14) 

1 - 0.47 [1 - exp(-1.5051z 0 - 4789 )] , for z t 0.4 


where z is a physical parameter defined by 

« . 4«, „ 4^ „ 3 (DM (15) 

r\ 4 (D/V) ' 


The coefficients in Eq . (14) are: f. - 1.1128, f 2 
- -2.5576, f 3 - 5.6257, and f 4 - -5.9498. ‘ 


4. APPLICATION TO SATELLITE REMOTE SENSING 

Using these parameterizations , the scatter- 
ing and absorption properties for ice crystal size 
distributions obtained during the FIRE Intensive 
Field Observation and the effects of small ice 
crystals on the infrared radiative properties are 
examined. As shown in Fig. 1, the absorption 
efficiencies of ice crystals for the wavelengths of 
11 and 12 pm can differ by as much as a factor of 
two. This could provide a means for the detection 
of small ice crystals using the radiances measured 
from AVHRR channels 4 and 5. 
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Fig. 4: Single scattering albedo for randomly 
oriented (a) oblate and (b) prolate spheroids with 
the axial ratio of 2 at a wavelength of 11 pm 
computed from the large particle parameterization 
(solid line), the small particle parameterization 
(dash-dotted line) , and the exact method (circles) . 





